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INTRODUCTION 

This  report  presents  the  mathematical  development  of  coupled  equations 
whose  solutions  are  expected  to  give  the  variation  in  pressure  along  the 
axis  of  tip  vortices  commonly  generated  by  lifting  surfaces  (hydrofoils, 
propeller  blades,  etc.)*  This  work  was  motivated  by  the  long-observed 
lack  of  agreement  of  cavitation  inception  speeds  measured  with  models  in 
water  tunnels  and  those  determined  at  ship  scale.  Generally,  the  model 
prediction  is  highly  non-conservative,  yielding  cavitation  inception  speeds 
which  are  roughly  twice  those  observed  during  prototype  operation. 

This  lack  of  scaling  of  tip  vortex  inception  speed  is  countered 
through  the  use  of  empirical  correction  formulae  which  are  based  on  the 
195**  work  of  McCormick1  who  resorted  to  a  semi-empirical  digest  of  data  which 
was  ultimately  published  in  1 9622 .  Somewhat  more  recent  (1961)  theoretical 

q 

work  is  that  of  M.G.  Hall  who  developed  a  theory  stemming  from  fundamentals 
which  shows  quite  good  correlationwi th  pressure  distributions  in  delta  wing 
tip  vortices  in  air. 

A  highly  significant  outcome  of  Hall's  analysis  is  the  large  contribu¬ 
tion  to  the  pressure  reduction  along  the  vortex  axis  arising  from  the  in¬ 
duced  axial  velocity  component.  McCormick's  empiricism  does  not  include 
this  component.  However,  significant  differences  between  Hall's  theory 
and  measurements  remain.  The  current  analysis  is  directed  at  overcoming  a 
deficiency  of  Hall's  theory  which  employed  an  outer  solution  in  which  the 
velocity  and  pressure  are  only  functions  of  the  angular  variable  and  inde¬ 
pendent  of  the  radial  coordinate.  As  the  solution  must  tend  to  uniform 
flow  for  large  radial  distances,  it  is  not  possible  to  reconcile  Hall's 
assumption.  Hence,  the  need  for  a  more  general  formulation  seems  apparent. 

The  theoretical  model  developed  herein  regards  the  flow  as  viscous 
everywhere  in  contrast  to  many  previous  studies  which  assume  at  the  outset 
that  the  flow  is  potential  except  for  a  viscous  core  region.  Inner  and 
outer  flow  regimes  are  encountered  giving  rise  to  a  matching  procedure 
along  a  conical  surface  having  a  very  small  apex  angle.  The  swirling  motion 
is  generated  by  an  input  transverse  velocity  which  is  obtainable  from  the 
solution  for  the  flow  about  the  lifting  surface  of  interest.  The  resulting 
coupled  equations  are  complicated,  requiring  extensive  computational  efforts 
to  produce  numerical  results. 
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This  study  was  supported  by  Contract  N00014-78-C-01 14  under  the 
Nava)  Sea  Systems  Command  General  Hydromechanics  Research  Program,  tech¬ 
nically  administered  by  the  David  Taylor  Naval  Ship  Research  and  Develop¬ 
ment  Center,  Bethesda,  Maryland.  Davidson  Laboratory  designation  was 
Project  040. 

1.  Statement  of  Problem 

We  assume  that  there  exists  a  conical  vortex  of  extremely  small  apex 
angle  whose  vertex  Is  at  0  and  whose  axis  of  symmetry  is  along  the  z-direction, 
which  is  also  the  direction  of  a  uniform  stream  at  infinity.  We  take  cylin¬ 
drical  coordinates  (r,0,z)  and  the  complete  velocity  field  will  be 

ur  +  v0  +  (W+w) z  ,  (1.1) 


W  being  the  uniform  stream  at  infinity.  The  flow  is  axially  symmetric  and, 
on  the  basis  of  "Oseen's  approximation",  the  equations  of  motion  and  con¬ 
tinuity  will  be  taken  in  the  following  form: 


W  — 
3z  r 


v(y2u  _  4) 
p  3r  rz 


wll 

3z 


v(v2v  -  -j^r) 


W  ~ 

3z 


=  -  1|£  +  v  (vzw) 
p  a  z 


J7  (ru)  +  Ti  (rw)  =  0 


where 


V2  3  -li  +  + 

3z  3r2  r  3r 


(1.2) 

(1.3) 

(1.4) 

(1.5) 

(1.6) 


is  the  Laplacian  and  v  is  the  kinematic  viscosity. 

Let  n  be  the  vorticity  about  the  0  direction,  so  that 


3u  3w 
n  "  Fz  '  3? 

If  we  eliminate  the  pressure  p  between  (1.2)  and  (1.4),  we  obtain 


(1.7) 
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"t'  ''HI  ’  i  tf-r* 


WU -|vvz.v(^-^) 


From  (1.5)  we  can  introduce  the  stream  function  i|i  such  that 


3il>  3i^ 

ru  ■  »  .  rw  *  —  -  ■ ' 

3z  ’  3r 


and,  from  (1.7),  we  have 


rn  ■  4i  _  +  \li  -  —  ill 
zz  rrr  r  yr 


(1.10) 


If,  in  place  of  i|>  ,  we  use  the  function  V  defined  by 


ifi  -  rY 


(1.11) 


n  « 


(1.12) 


We  note  now  that  Equation  (1.3)  is  a  p.d.e.  for  v,  Equation  (1.8)  connects 
n  and  v.  Equation  (1.12)  connects  the  function  V  (and,  therefore,  the 
stream  function  iji)  with  n,  and  the  problem  consists  in  solving  these  equa¬ 
tions  successively. 

We  are  assuming  that  there  is  no  vorticity  in  the  oncoming  stream 
Wz  and  in  solving  (1.8)  we  shall  impose  the  condition 


n  =  0  ,  z  =  0 


(1.13) 


and  we  shall  solve  for  n  in  z  >  0.  Likewise,  in  solving  (1.12)  for  y, 
we  shall  impose  the  condition 


y  =  0 


z  =  0 


(1.14) 


and  we  shall  solve  for  t  In  z  >  0.  The  conditions  which  will  be  imposed 

upon  v  will  be  discussed  later,  but  we  envisage  that  the  vortex  is  driven 

by  an  imposed  transverse  inflow  which  is  known  along  a  generator  of  the 

cone  r  *■  z  tan  0  . 

o 

2.  Inner  and  Outer  Solutions  for  the  Transverse  Velocity  v 

It  Is  convenient  to  introduce  the  parameter  k  defined  by 

2k  *  — 

M 


> 


3 


(2.1) 
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and  (1.3)  can  then  be  written  in  the  form 


72v  -  —jr  ~  2k  =  0 
r^  3z 


If  we  write 


then,  using  (1.6) 

93X  +  33X  +  _1_  32X 

3zz3r  3r3  r  3rz 


1  3X 

r7  3r 


2k  =  0 

3z3r 


that  is 


J_  +  lii  +  I  ax  .  2k  ax,  u  n 

3r  3z^  Jr7  r  3r  3z; 


and  we  take  X  to  satisfy  the  equation 


ifx+3fx+I3X_2k3Xo0 
3z2  3r2  r  3r  K  3z 


If  we  write 


kz 

X  =  e  4> 


then  $  will  satisfy  the  p.d.e. 


>4  ♦  14  ♦  i  !* .  „*♦ .  o 

dzz  dr^  r  3r 


(V2  -  k2)$  =  0 


(2.6) 


The  functions  v  and  $  will  be  related  as  follows: 


kz  3$ 


We  can  now  build  up  solutions  for  0  using  the  method  of  separation  of 
variables.  If  we  look  for  a  solution  for  <J>  of  the  form 


<J>  =  e  mZf(r) 


k 
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then 


(m2  -  k2)f  =  0 


(2.9) 


There  are  two  cases  to  discuss  according  as  m  >  k. 

Case  1  m  >  k:  Inner  Solution 

In  this  case,  the  acceptable  solution  of  (2.9)  which  is  finite  at 
r  =  0  is 


f(r)  =  Jo{rt42-k2}  (2.10) 

and  thus  we  can  build  up  a  general  solution  for  4>  of  the  form 

oo 

<f.(r,z)  =  -  /  p(m)e  mzJ  { r/m2  -  k2 }dm  (2.11) 

i  ,  o 

m=k 

where  p(m)  is  an  arbitrary  function  of  m. 

Case  2  0  <  m  <  k:  Outer  solution 

In  this  case,  we  use  the  solution  of 

-  <k2-”2)f  -  0  (2-12) 

which  tends  to  zero  as  r  namely 

f(r)  =  K  { r/k2  -  mz }  (2.13) 

o 

We  can  now  build  up  the  general  solution  for  4>  in  the  form 
k  _  _______ 

4>e(r,z)  =  -  /  v(m)e  mZ  K  {r/kz  -  m2}dm  (2.H) 

m=o  ° 

The  formulae  for  v  now  follow  from  (2.7)  and,  using  J'(ii')  =  -  J,(i>), 

o  * 

10  (iji)  =  -  we  obtain 

OO 

v.(r,z)  =  e^2  /  p (m) /m2  -  k2  e  mZ  J  { r/mz  -  k2 )dm  (2.15) 

'  m=k  1 


R-2 1 24 
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m  =  k  cos0 

(2.2k) 

B{0)  =  k2sin20  v(k  cos0) 

(2.25) 

then  the  integral  equation  for  B(o)  will  be 

/  B(0)  e"kR  COs9  COs9o  K  { kR  sin0  sin6  }d9  =  V  (R)e"kR  COS0o  (2.26) 

1  loo 

o 

It  is  clear  that,  when  we  proceed  to  calculate  v.(r,z)  and  vg(r,z)  in  the 

above  way,  these  two  solutions  will  be  continuous  on  the  cone  6  *  0q; 

however,  the  normal  derivatives  3v./3n  and  3v  /3n  will  be  discontinuous  on 

i  e 

0  =  0  •  The  sketch  below  defines  the  role  of  the  driving  transverse  velocity 
produced  by  the  lifting  surface  along  the  generator  of  the  matching  cone  from 
the  leading  edge  to  downstream  infinity. 


I 

;  * 

i 
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3.  The  Solutions  for  the  Vortioity  n  and  the  Stream  Function  tp 

We  can  write  (1.8)  in  the  form 

v2n  -  -4-  -  2k  f2-  =  -  —  vv 

rz  3z  \>r  z 


(3-D 


If  in  (3.1)  we  wri te 
kz  3a 


n  *  e 


we  obtain 


3r 


(3-2) 


—  (a  +  a  +  —  a  k2a)  = 
3r  zz  rr  r  r 


2  -kz 
—  e  vv 
vr  z 


(3-3) 


Thus  we  can  take  a  to  satisfy  the  equation 


(V2  -  k2)a  =  f  (r,z)  (3.14) 

where 

fM-.i/I  .  (3.5) 

o 

We  can  always  add  to  a  an  arbitrary  function  of  z,  say  A(z),  but  this  will 
leave  n  unchanged  and  thus  we  can  take  A(z)  =  0.  The  principal  requirement 
upon  n  is  that  it  should  tend  to  zero  at  infinity  and  thus  we  can  look  for 
the  solution  for  a  which  satisfies  (3  -  ^ )  and  which  tends  to  zero  at  in¬ 
finity.  The  solution  of  this  problem  is  obtained  in  the  Appendix  and  is 
given  by 

00  °°  2-ir  -kR 

i4Tra(ro,z  )  =  -  /  |  /  |-  f(r,z)r  dr  dz  d$  (3.6) 

0  z=— 00  r=o  <f>«o 

where 

R2  =  r2  +  r2  -  2rr  cos<Ji  +  (z-z  )2  (3-7) 

00  o 

It  follows  from  (3-2)  that  the  function  n  is  given  by 

kz  03  00  2rr  -kR 

i»irn(ro,zo)  =  -e  0  f  f  J  f  f(r,z)r  dr  d<p  dz  (3-8) 

0-0  o  o 

We  can  simplify  the  r.h.s.  of  (3.8)  as  follows;  we  write 
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CM 

8 

-kR 

J  = 

/  / 1 
o  o 

f(r 

we  have 

3J 

3r 

CO  2lT 

„  -I 

"/  / 

3  /  e 

Sr-  R" 

o 

o  o 

o 

*= 

CM 

8 

.  -1 

-/  / 
o  o 

3  ,e 
3R  R 

00 

2ir  ^ 

-kR 

-/ 

1  — 

1  dR 

(-  )  ■ 
'R  ' 

o 

o 

CM 

8 

-/  / 
o  o 

d  ,e 
dR  VR 

kR  (r  -  r  cos<f>) 


r  f  (r,z)dr  d<{> 


rf(r,z)dr  d(J> 


3R  s  >  n<[>  3R, 
Jr  r  3^ 


(3.9) 


Hence 


3  J 
dr 


<=°  2n  -kR  .  .  -kR 

=  /  /  {‘cos^(t  >  +-T±'k  (f — )}  r  f(r’z)dr  d* 

o  o  o  Y 


(3.10) 


We  can  now  integrate  by  parts  the  terms  on  the  r.h.s.  of  (3-10)  to  give 


„  2tt  -kR®  00  2tr  -kR  . 

=  -  /  [r  coS(f>  f  (r,z)  f-  ]  d<f>  +  /  J  cos  —  (rf  (r  ,z)  }dr  d£ 


o  o 


-kR  0  0  -kR 

00  2tt  a>  2tt 

+  /  [s in(j>  f(r,z)  ]  dr-/  /  ^  f(r,z)  —  (s i n<t>)dr  d(j> 

o  4>=o  o  o 


Since  there  is  no  contribution  from  the  square  bracket  terms,  we  have 

co  2ir 


_  J- D 

■—  =  /  /  cos<|)  (-~(rf)  "  f)dr  d<fl 

o  o  o 


7  |2tT  ,  e  3f(r,z)  rdrdiji 

/  J  cos*  ¥ 

o  o 


(3.11) 
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It  then  follows  from  (3.11)  and  (3-8)  that 


,  ,  _  x  kZo  7  "  2?  -kR  .  3f(r,z)  ,  , 

*»TTn(ro,zo)  =  -e  J  /  j  e  cost  - ^7— ^  r  dr  dz  d<j> 

-»  o  o  R 


(3.12) 


and  we  can  now  substitute  for  3f/3r  using  (3.5),  namely 


3f  2  -kz  /  x 

-=--e  v(r,z) 


(3.13) 


so  that  we  have 

kz  00  °°  2tt  -kz-kR  ,  . 

27rvn(ro>zo)  =  e  °  /  /  /  - — ^ -  cosi{>  v(r,z)  -■->*!■  dr  dz  d?  ( 3  - 1 

and  we  note  that  the  formula  (3*14)  relates  n  directly  with  the  transverse 
velocity  v. 

We  consider  next  the  formula  for  the  stream  function  in  terms  of  n, 
namely  the  equation  (1.12).  If  in  (1.12)  we  write 


(3-15) 


we  obtain 


F  ‘  6rr  +  7  6r»  '  "<r'z> 


(3- 16) 


so  that 


V26  =  F(r,z) 


(3.17) 


where 


F (r , z)  =  /  n(o,z)do 


(3.18) 


The  solution  for  0  which  vanishes  at  infinity  is  given  by 


CO  00  2tT 

|,7rB(ro,zo)  =  -  /  /  / 

-00  o  o 


F(r  ,z)  r  dr  d<|>  dz 
R 


(3-19) 


where  R  is  defined  in  (3.7).  It  follows  from  (3 . 19)  and  (3.15)  that 


R-21 24 


i*7t4' (r0*z0)  =  "  aF"  /  /  /  F  F(r,z)r  dr  d<J>  dz 
o  -«°  o  o 


(3.20) 


If  we  consider 


I  —  /  /  p-  F(r, z)r  dr  d<j> 


(3.21) 


we  have 


JT=  /  /  sHlf*  F(r»z)r  dr  d<t> 

o  o  o  o 


00  2ir  (r  cos<(!  -  r  ) 


F(r  ,z)r  dr  d<f> 


00  2ir  cos4>(r  -  r  cos<f>)  -  r  sin24> 

/  /  - ^“51 - 2 -  F  (r  ,z)  r  dr  d<j> 


OO 

=  /  /{-COS4.  ~(£)  +  F(r,z)  rdrd^ 


o  o  o 


(3.22) 


We  integrate  by  parts  the  two  portions  of  the  integral  (3-22)  and  we  obtain 


=  7  [.  £22±  r  F(r  .*)]""%♦  ♦  /  7  (r  F(r,z)}dr  d* 


o  o 


r=o  o  o 


./(ii£i.rF(r,2))2'  <r.  » (si„*)drd4 

o  4>=o  o  o 

Provided  F(r,z)  tends  to  zero  at  r  -*■  «,  which  will  be  satisified  in  the 
present  case,  there  is  no  contribution  from  the  square  brackets;  hence 
we  obtain 


“  «  J  /  ^  (£(rF)  -  F>dr  <*<*■ 

ro  o  o 


-/  / 


r  cosj>  5F(r ,z) 


dr  d<j> 


(3.23) 


» 
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/  /  ^  n(r ,z)cos*  r  dr  d* 


(4.5) 


31 1  "  2it  a  i 

tt-  =  /  /  tt"  (4-)  n(r  ,z)cos*  r  dr  d<j> 


9r  1  1  3r  'R 

o  o  o  o 


00  2tt  r  cos*  -  r 


n  cos*  r  dr  d* 


00  2ir  cos*(r  -  r  cos*)  -  r  sin2* 


ncos*  r  dr  d* 


*  /  /  {-cos<t>  •^•(■jjr)  n  cos*  r  dr  d* 


/  /  {-  cos2*  +  sin*  cos *  n(r,z)  rdrd<()  (4.6) 


We  can  now  integrate  by  parts  the  two  portions  of  the  integral  in 
(4.6)  to  give 


3 1  2tt  °°  °°  2it  2*  -n 

~  I  ("rr'  cos2*  •  ^]  d*  +  J  /  C°~R  Jf  (rn)  dr  d* 
o  o  r=o  o  o 


w  1  411  3 

+  /  [n  sin*  cos*  •  dr  -  /  /  —  (sin*  cos*)dr  d* 

o  *  *  =o  o  o  * 


and,  since  there  is  no  contribution  from  the  square  brackets,  we  have 


o  I  w  L\\  2x  A 

r~  -  /  /  ^  tt  (rn)  -  £  (cos2*  -  s  in2*)  }dr  d* 


3r  >  1  R  3r 

o  o  o 


=  /  /  4-  (cos2*(r  ~  +  n)  -  n(cos2*  -sin2*)}drd* 
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/  /  jj-  {r  cos2<j>  +  n  sin2$}drd$ 


(4.7) 


It  now  follows  from  (I*. 4)  and  (4.7)  that 


3H* ( r  ,z  )  ®  ®  2n 


"  /  /  /  r-  (r  cos2$  37  +  n  sin2$}dr  dzd<|i  (4.8) 

-•  o  o  0 


this  result  being  valid  for  ail  r  and  z  . 

o  o 

We  can  deduce  from  (4.8)  the  limiting  value  of  3H'/3ro  as  r  0;  using 
(4.1),  we  have 


4n  11m  -  °  •••° 

r  -*>  3ro 

o 


3f(r  ,z  )  •»  ®  2n 

1  J  f 


r  cos 2<ft  3n/3r +  n  sln2d 


ir  dz  dt|> 


'r  +  U-z 


Hence 


3f(r  ,z) 


4it  1  im  - - 

r  -o  3ro 

o 


r  3n/3r  +  n 


dr  dz 


-«  o  rr£ 


(4.9) 


Using  (4.2),  we  also  have  the  result 


3'f (r  ,z  )  >f(r  z  ) 

lim  - r - *  lim  - 

dr  r 

r  -k>  o  r  -*o  o 

o  o 


(4.10) 


so  that  from  (4.9)  we  deduce  that 


,  |m  {.lllo.V.  +  y(r°,-Z-°.)l}  *-1-77  9/ij(rn)drd_2 

r  -*o  ro  ro  o  /r2  +  (z-z  )2 

o  o 


(4.11) 


Using  the  stream  function  'i'  and  the  results  (1.9),  (1.11),  we  have  the 

general  formula  for  the  velocity  component  w(r  ,z  )  in  the  z-direction, 

o  o 

3f(r  ,z  )  »(r  ,z  ) 

w(r  ,z  ) - r?-2 - 2_2_ 


From  (4.11)  and  (4.12),  we  obtain  the  result 


(4.12) 
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w(0.xo)  -  Urn  w(ro.ro)  -  J.  /  / 


ro’° 


-  o  Vr*  +  (z-z  )z 
o 


(4.13) 


We  note  from  (4.13)  that  w(0,z  )  vanishes  as  z  -*•  +  •».  We  can  obtain  an 

o  o 

alternative  formula  for  w(0,zQ)  if  we  use  integration  by  parts  in  (4.13). 
This  gives 


w(0,z  )  m  J  /  I"  -rri^Z^  ■  ]  dz  -  j  }  I  rr )  •£-  {■  ■  —  1  —  ■  }drdz 

/rz  +  (z-z  )z  r*o  o  /r2  +  (z-z  )z 


Hence 


w(0, z0}  "  I  /  / 


r2n(r ,z)drdz 

o  {r2  + (z-z  )2)3/2 
o 


(4.14) 


and  this  result  makes  it  clear  that  w(0,z  ) ‘vA/z  3  when  z  is  large. 

o  o  o 

By  making  use  of  (3*14)  and  (4.14),  we  can  express  w(0,zq)  in  terms 
of  the  transverse  velocity  v.  With  a  slight  change  of  notation,  we  have 


+oo  »  r  2  n(r  ,z  )dr  dz 

w(o,?)  -  i  /  /  — — s--a-  -2-f 

V*  rt,'°  {ro2  *  <V£>2' 


so  that 


kz 


4irv  w(0,£)  -  /  / - — yjj 

*.— •  r  -o  (r  2  +  (z  -£)zr  / 


r  2  e  °dr  dz  «  ®  2ir  -kz-kR  „  ,  . 

- — — °  /  /  /  - — = - cosi}>v  (r  ,z)  !~-dr  dg> 


o  o  o  o 


-oo  o  O 


3z 


and  by  interchanging  orders  of  integration,  we  can  write  this  in  the  form 


OO  OO  2lT  /  v 

4irv  w(0,C)  *  /  /  /  e  Zcosif>v(r  ,z)  3  3^*—  «(r  ,z,4>  ,£)dr  dz  d<f>  (4.15) 
-00  o  o 


where 


o>(r,z,4i,C)  *  /  / 


kz  -kR 

*  r  e  dr  dz 

o  00 


i  i  - 579  * 

z  *-«  r  -o  R{r  2  +  (z  -C)2r 


(4.16) 


In  (4.16)  R  is  defined  by  (3-7),  namely 


R2  *  r2  +  r  2  -  2r  r  cos$  +  ^z-z  )2 
00  o 


(4.17) 
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5.  Tiie  Determination  of  the  Pressure  Function 

It  is  easily  shown  that  the  equations  (1.2)  and  (1.4)  can  be  written 
in  the  form 


I  „  w  2i L.  I2  .21  (rn) 

p  3r  3z  r  r  3z  'rn; 


(5.1) 


;  If 


(5.2) 


Elimination  of  p  between  (5.1)  and  (5.2)  gives  rise  to  the  vorticity  equation 
(1.8);  accordingly  the  right  hand  side  of  the  equation 


“  7dp  =  “  7(lf  dr+!f  dz)  "  {wfl-F2_7^(rn)}dr  +  {wl7+7^  (rn)}dz 


(5.3) 

is  necessarily  a  perfect  differential  and,  in  theory,  we  should  be  able 
to  determine  p  from  the  line  integral 

'f  ’  {  {“lf-7  ‘7  &'">>«  *  <w  I?  *  7  fr'Mo*  «-4> 

As  an  alternative  approach,  we  note  that  (5.1)  and  (5.2)  can  be  written 
in  the  form 


7  If  =  w£(ru)  •  v2  '  v£  (rn) 


(5.5) 


‘  7  If  =  W£(rw)  +  v  TF  (ril) 

Differentiating  (5-5)  w.r.t.  r,  (5.6)  w.r.t.  z  and  adding,  we  obtain 

-  7  (i<r  !?>  ♦• r  ■“  n-  <£<™> +  >  -  2“  If  • 

and,  using  (1.5)  ,  this  becomes 

L  {i!£  +  1  i£  +  ill}  =  2v  li 
p  '■Jr7  r  3r  Tz71  3z 


(5.6) 


or 


1  _2  2  3v 

7  7  p  ”  7  v  37 


(5.7) 
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A  particular  solution  for  p  of  (5-7)  is  given  by 


"  ■  2n 


p"{ro'zo)m"fa  !  f  / 


2/r  v(r,z)  3v(r  ,z)/3z  rdr  dz  d4 


If  we  write 


z*-<®  r*o  o  /rz+rz~  2rr  cos4>  +  (z-z  ) 
o  o  o 


p  *  p  +  P 


then  P  will  be  a  harmonic  function,  i.e.,  satisfying 


(5.8) 


V2P  *  0 


(5.9) 


The  pressure  function  p  is  continuous  throughout  the  space  occupied  by  the 
liquid,  but  at  the  surface  of  the  cone  0  * 0q  (see  Section  2)  there  will  be 
a  discontinuity  in  3p/3n  where  n  is  the  normal  to  the  surface  of  the  cone. 

fc  Accordingly,  in  order  to  deter- 

✓  >/  \  mine  P,  we  must  consider  an  in- 

x  ^  ""  f.  '  ~ 

\  ner  solution  P.  valid  in  O<0<0 

- L-V,  4. 


and  also  an  outer  solution  P 

— € 

valid  in  0  <  0  <  it  ,  with 
o 


P.  =  P 
—i  — e 


0-0 


(5.11) 


a  4  0  at  0  =  0 


(5. 12) 


This  pattern  of  behavior  of  P^  can  be  achieved  by  distributing  pressure 

sources  over  the  cone  0=0.  The  function  a  entering  into  (5-12)  will  be 

o 

determined  from  the  discontinuity  in  3v/3n  at  0  *  0Q,  v  being  the  transverse 
velocity.  Thus,  the  determination  of  p  using  (5-7)  will  lead  eventually 
to  an  integral  equation  for  the  unknown  pressure  source  distribution  on  0-e^. 

If  we  are  interested  solely  In  the  pressure  distribution  along  the 

vortex  axis  r  =0,  z  >0,  it  is  possible  to  avoid  the  solution  of  an  in- 
o  o 

tegrai  equation  by  returning  to  (5-3).  If  we  take  rQ  -0  and  dro  =  0  in 
(5- 3) ,  we  obtain 


.  3w(r  , z  )  3n(r  ,z  )  n(r  ,z  ) 

- dp  ■  (W  - — -  +  v  ( - — -  +  — — - )  1  dz 

p  3z  3r  3r  o 

o  o  o  r  *o 

o 


(5.13) 
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Thus  p(0,z  )  will  be  given  by 


p(0,z  J  p, 


+  -  =  W  w(0,z  )  +  v  /  {- 


zo  3n(r  ,z)  n(r  z  ) 


3r 


r  -}  dzQ  (5.1*0 

o  r  =o 
o 


In  order  to  determine  the  integrand  of  the  integral  in  (5.1*0.  we  return 
to  Eq.  (3-8)  from  which  we  obtain 


4*  +~)  =  ~e  °  (g~z-+ j-  ^_)  /  /  /  f  f  (r,z)r  dr  d<f>  dz  ,  (5.15) 


kzo  ,  3*  1  3  -  -  2*  _-kR 


o  o  -°°  o  o 


f(r,z)  being  defined  in  (3.5).  from  the  Appendix,  we  see  that  if 


00  00  2tt  -kR 

^nr$(ro,zo)  =-/  /  /  |-  f  (r  ,z)  r  dr  d$  dz 

-00  O  O 


(5.16) 


then 


Ji*  +_L_iL  +  Ji*  k2*  =  f(r  z) 

JT7  r  3r  TT7  1  o*V 

o  o  o  o 


(5.17) 


If  we  apply  this  result  to  the  right-hand  side  of  (5.15) »  we  obtain 

kz  »2  oo  oo  2ir  -kR  kz 

l*ir  +  =  +  e  °  (gj-r  *  k2)  f  f  f  ^  f  (r  ,z)r  dr  d$  dz  +  Aire  °f(ro,zo) 

O  o  O  -00  o  o 

(5.18) 

It  is  permissible  to  make  rQ  -*•  0  in  the  result  (5-18)  and,  if  we  make 
use  of  the  definition  of  f(r,z)  in  (3.5),  we  obtain 


wg-'r* 

o  o 


kz 


o,  3Z 


oo  co  2tT 


(jfr-K2)/  /  / 


-k*Vz+(z-zo)z 


r  =o 
o 


o  o  /rz  +  (z-z  )^ 


f  (r  ,z)r  dr  d$  dz 


kz 


-  2.e  °  (#t-k2)  /  /  e 


-k/rz  +  (z-z  )2 


3z 


-oo  o  /rz  +  (z-z  )z 
o 


f  (r  ,z)  r  dr  dz 


(5.19) 


Thus  (5*  I1*)  can  be  written  in  the  form 
p(0,z  )  p 


zo  kz  „  -kA‘:  +  (z-z0)z 

_2_+-J!  =  W  w(0,zo)  +jv  J  e  °f(~7--k2)  /  /  -  f(r,z)r  drdz)dzo 


o  /r^  +  (z-z  ) ^ 
o 


(5.20) 
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and  w(0,zq)  is  given  by  (4.14)  or  (4.15).  Since  we  have 

Zo  kz  a2  kz  j 

/  e  -  k2a)dz  =  e  °  (^_  -  k«)  , 

“  o  o 

provided  a(zQ)  and  a* (zq)  tend  to  zero  at  infinity,  it  follows  from  (5.20) 
that 


p(0,z  )  p.  .  kz„  -  »  -k<'r^(2-zo)^ 

- - - *—  .  H„(0,zo)+jve  (jf--k)/  i  —  -  —  f(r.z)rdrdz 

o  -o  /r<i*(z-zo)“  (5.21) 

We  have 


-k/r  + (z-z  )2 


re 


-  f(r,z)dr  =  -  —  [  —  { 

'  r~7 1 t - nr  k  J  3r 

o  /r2  +  (z-z  )z  o 


,  “  -k/r2+(z-z  )2 

1  r  3  ,  o 


}f (r ,z)dr 


i  -k^*(z-z0i«  f(r.z)]“  ^  ,  J^lz-z^ 

k  1  r=o  o  3 


dr 


and,  using  the  definition  of  f(r,z)  in  (3-5),  we  see  that  the  square  bracket 
makes  no  contribution  and  thus  we  have 


-k/r2  +  (z-z  )z 


re 


°  /r2+(z-z  )2 

o 


oo 

f (r,z)dr  =  r-  / 


-k/r2+(z-z  )2 


9'f~aV,Z'' dr  (5‘22) 


From  (3.5).  we  have 


=  .  i  i  e'kz  v (r  z)  3v(r>z) 
3r  v  r  vv  ,z  1  3z 


(5.23) 


hence  (5-21)  can  be  written  in  the  form 


P(0,Z  )  p  ,  kz  .  oo  oo 

- +  —  =  W  w(0 , z  )  -r-  e  ° (-5 - k)  /  /  —  e 

P  p  o  k  vdzo  •  Jq  r 


-k/rz+(z-z  )z 
o 


■k2.(r,z)^%^drdz 

(5.24) 


The  result  (5.24)  is  valid  along  the  vortex  axis  r  =  0,  z  >0.  A  general 

o  o 

equation  for  the  pressure  p  can  be  derived  from  (5.4)  as  follows.  We  write 
(5-4)  in  the  form 
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-  £  ’  /  <W  If  r  ■vlt)dr  *  tU  H  *u  <1?  *  f),d2 


3w  A  ..  ,3n  ^  n\ 


-  /  CW(2J  -  §*)dr  -  f  dr  -  vfj  dr  +  W(§*  dr  +g  dz)  +v  (f?  ♦  ?)dz) 


dw  .  , 3w 


rln  ^  ru 


Hence 


-  ■£■  =Ww(r,z)  +  /  {w  ndr  -  j  dr  -  v~  dr  +  v  (f7  +  ~)dz) 


(5-25) 


It  is  convenient  now  to  write  n  =>  e  z  3a/3r  as  in  (3-2),  where  the 
function  a  satisfies  (3-1*)-  We  then  have 


-  f  .  Ww(r,z)  ♦  /  (W  ekk  £  dr -fdr  -  «  £  (.»  £>dr  ♦«"(£♦*)  £  dz) 

^  r' 


3  /  kz  3av 


kz,  3  ,  1 \  3a 


and  since,  by  the  definition  (2.1),  W  =  2kv,  this  becomes 


-  £  .  W„(r,z)  +  /(2kw  ak*  |£  dr  -f  dr  -v  e“«jgL*k  =»dr  £  dz) 


kz,  32a  ,  .  9a> 


kz,  3  ,  1 .  3a 


=  Ww(r  ,z)  +/{vekz(~7-  +  7  '  vekz(~7  -  k  |y)dr  -  £  dr}  (5.26) 


The  function  a(r,z)  satisfies  the  differential  equation  (3-^),  namely 


l!«  +  1  |SL  +  1^.  -  k2a  =  f  (r  ,z) 
3r2  r  3r  3zz 


(5.27) 


and  we  replace  the  coefficient  of  dz  in  (5-26)  using  (5.27)  to  give 


32a  ,  ,  ,  .  kz  3  /3a 


-£=Ww(r,z)  +  /  (e  Z[f(r,z)  ~jpr  +  k2a]dz  -  ve  77(77-  ka)dr  --  dr} 


(5.28) 


We  have 


±  -  to)  -  ek‘(f^  -  A) 


Hence 
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-  £  =  Ww(r,z)  +  /  {v  ekZf  (r,z)dz  -y  dr  -  v -^-[ekz  (|j  -  ka)]dz  -  ve  ar  ^3z  '  ko^dr* 


=  Ww  (r,z)+/  (ve  Zf  (r ,  z)dz  -  p- dr  -  vd  [e  Z(— --ka)]} 


=  W  w(r  ,z)  -  vekZ  (|—  -  ka)  +  /{veKZ  f(r,z)dr  -p-  dr) 


We  now  introduce  a  function  y(r,z)  having  the  property 


ekzf(r,z)=|l 


(5.29) 


(5-30) 


where  f  is  defined  in  (3-5) 
r ,  we  obtain 


kz  3_f_  _  32y  . 

e  3r  3r3z 


If  we  differentiate  (5-30)  with  respect  to 


replacing  3f/3r  from  (5.23),  we  have 


Jh-  .  .  1  1  v(r  z)  1V-(I’Z) 

3r3z  v  r  V'  ’  '  3z 


Hence 


32y  =  _  J_  , v2 s 
drdz  v  3z  V 


and  we  take 


3y  J_  v* 

3r  v  r 


(5.31) 


If  we  use  (5.30)  and  (5-31)  in  (5-29),  we  obtain 

-  ■£•  =  Ww(r,z)  -  vekz  (■—  -  ka)  +  /v(-|£  dz  +  dr) 


giving 


-  £  =  C  +Ww(r,z)  -  vekz(|^-  -  ka)  +  vy(r,z) 

p  O  o  Z 


(5.32) 


where  C  is  a  constant  and  the  function  y(r,z)  will  be  given  by  the  line 
o 

integral 
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Y (r , z)  =  /  {ekzf(r,z)  dz  -  dr} 
c 


(5.33) 


If  we  choose  the  initial  point  of  the  curve  C  to  be  at  infinity 


so  that 


P(r.z) 


y(r,z)  *  /  (ekz  f(r,z)dz  -  dr} 


(5-34) 


then  the  functions  w,  eKZ(a  -ka)  and  y  will  tend  to  zero  at  infinity  and 


the  constant  Cq  will  be  given  by 


C  =  "P  /p 
o  00 


(5-35) 


It  will  be  noted  that,  when  the  pressure  along  the  vortex  axis  is  deter¬ 
mined  from  the  general  formula  (5-32),  we  have  r  =  0  and  dr  =  0  in  (5.3*0 
and,  since  f  =  0  on  r  =  0,  it  follows  that  Y(r>2)  =  0;  the  formula  (5-32) 
then  gives  the  result  (5.24).  The  result  (5.24)  is  the  only  simple  formula 
for  the  pressure;  in  all  other  cases  it  is  necessary  to  determine  the  func¬ 
tion  y(r,z)  defined  in  (5-34).  When  we  substitute  in  (5-32)  for  the  func¬ 
tion  a(r,z)  defined  in  (3-6),  we  obtain  the  following  general  formula  for 
the  pressure 


-p  o  .  +“>  “>  2tt  -kR 

"r  =  Ww(ro'z0)+¥  (TT"k)  /  /  /  f  f(r,z)r  dr  dz  d  *  +  vy(r,z) 

O  z=-oo  o  o 


(5-36) 
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SUMMARY 

We  summarize  the  many  results  presented  in  the  previous  sections 
by  pointing  out  the  numerical  work  required  in  order  to  have  an  under¬ 
standing  of  the  physics  underlying  the  mathematics.  The  same  numerical 
results  can  also,  of  course,  be  used  to  make  comparisons  with  pressure 
measurements  and  also  to  yield  predictions  with  regard  to  cavitation  in¬ 
ception. 

The  first  step  is  to  solve  the  two  integral  equations  for  A(0)  and 
B(6)  as  presented  in  Equations  (2.22)  and  (2.26).  In  order  to  be  in  a 
position  to  solve  these  equations,  it  is  necessary  to  specify  the  function 
Vq(R)  as  well  as  the  angle  6  .  These  two  inputs  must  come  from  an  outside 
source,  presumably  from  the  solution  of  the  inviscid  lifting  surface  theory 
near  the  wing  tip.  Nevertheless,  the  schemata  for  solving  the  equations 
can  be  specified  numerically  even  without  any  specific  inputs  available. 
Once  A(0)  and  B(0)  have  been  determined,  the  velocity  v(r,z)  can  be  cal¬ 
culated  using  Equations  (2.15)  and  (2.16) 

The  axial  flow  along  the  vortex  axis  can  then  be  evaluated  from 
Equation  ( ^4 . 15) >  while  the  pressure  along  the  vortex  axis  is  determined 
from  Equation  (5.24),  and  the  pressure  elsewhere  is  determined  from  Equa¬ 
tion  (5-36). 

It  would  be  of  some  interest  to  determine  the  circumstances  under 
which  the  pressure  is  minimum  along  the  vortex  axis,  for  it  is  the  minimum 
pressure  that  controls  cavitation  inception.  If  the  minimum  pressure  can 
be  shown  to  lie  along  the  axis,  then  the  relatively  simple  Equation  (5-24) 
would  control  cavitation  inception  rather  than  the  more  complicated  Equa¬ 
tion  (5-36). 
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APPENDIX 

The  solution  of  the  p.d.e.  (V2-k2)y  =  f(x>y,z) ,  with  f(x,y,z )  given  and 
with  the  condition  V  -*0  at  infinity 3  V2  being  the-  three-dimensional  Laplacian. 

The  function  'V(x.,y,z)  satisfies 

(V2-k2)4<  =  f  (x,y,z)  (1) 

^2  ^2 

where  V2  =  -r— 7-  +  7-  +  rr  and  y  0  at  Infinity. 

dX  dy  dZ 

Introduce  a  Green's  Function  G  with  with  properties 

(V2-k2)G  =  0  at  all  points  except  (x  y  z  )  (2) 

o  o  o 

G  =  ^  near  (xoY0zo) ,  R2=  (x-xq)2  +  (y-yQ)2  +  (z-zq)2  (3) 

G  ->  0  at  infinity  (I*) 

We  use  the  identi ty 


/// (XV2<p  -  <J,V2x)dT  =  //(X  !£•  "  <t>f“)dS  +//  (X  -  <t>  |^)do  ,  (5) 


outward  normal  from  x 

sphere  radius  center  (x  y  z  ) 


n  being  the  outward  normal, 
i.e.,  drawn  away  from  the  volume 
x.  We  identify 

X  =  ¥  ,  4>  =  G 

and  from  (1)  and  (2),  we  have, 
within  x 


XV2<j>  -  <t>V2X  =  ¥V2G  -  GV2f 

=  -  k2)G  -  G(V2  -  k2)f 

=  -  Gf 

Thus  the  identity  (5)  becomes 


//<*!£-  6  £>«  +  //<*£-  ®  -  -  ///Gf  dx 


ay 


3G 


3¥» 


3n 


3n 


3n 


3n 


(6) 


A- 1 


For  the  integral  over  a,  we  have,  using  (3), 


=  -Y  —  (— ) 

3RW 


1  3Y_ 
R  3R 


_  Y  \_  3Y_ 

R2  R  3R 


Hence 


lim  //  (Y  ||  -  G  |£)do  =  lim  //  (^-  +  7  |^)e2do>  -  4irY (xoyftzo)  (7) 


e-*o  a 


3n 


G-H3 


£  3R' 


o'o  o' 


From  (6)  and  (7)  we  have 


l,rf(xoyoZo)  +  //  f£  “  G  |£)dS  "  "///  Gf  dT 
S  T 

With  appropriate  safeguards  upon  the  behavior  of  Y  and  G  at  infinity,  the 
integral  over  S  will  tend  to  zero  as  S  moves  off  to  infinity  giving  the 
result 


IjttY (x  ,y  ,z  )  = 
000 


In  order  to  determine 


x  =  xQ  +  X,  y  * 


Then  Eq.  (2)  becomes 


G  +  G  +  G 
XX  YY  ZZ 


///  Gf  dt 

T 

G,  we  change  the  origin  to  (X0Y0ZQ)  bY  writing 

y  +  Y ,  z  =  z  +  Z 
'o  o 


-  k2G  =  0 


(8) 


(9) 


(10) 


If  we  look  for  the  spherically  symmetric  solution  of  Eq.  (10),  then  G 
will  satisfy 


32G  _2  _3G 
3R7  R  3R 


k2G  =  0 


01) 


where  R2  =  X2  +  Y2  +  Z2.  If,  in  (11),  we  write 


it  follows  that  g  satisfies 


(12) 


lL 


so  that 

/b\  a  kR  .  d 
g(R)  =  Ae  +  Be 

Since  we  require  G  +  0  as  R  •>  ®,  we  choose  A  =  0  and  then,  to  satisfy  the 
condition  (3) >  we  choose  B  =  1;  hence 


•-* 


Accordingly,  we  can  write  (8)  in  the  form 


(VoZo)  =  ‘///  R  f  dT 

T 


exp{ - k/(x-x  )"+(y-y^)^+ (z-z  ) 

4ir1f(x  y  z  )  =  -///  - - -  —  . . . -  °-  f  (x,y,z)dxdydz 


T  r(x-xQ;^  +  (V-yQr  +  (Z-ZQ) 


In  problems  with  axial  symmetry  where  x  =  r  cose,  y  =  r  sine,  xq  =  r^cose 

y  *  r  sine  ,  we  can  write  the  solution  for  4*  in  the  form 
7o  o  o’ 


^(rQ,zo) 


+“  «®  2it  -kR 

/  /  If  f  (r,z)  rdrdzd$ 


z=-°°  o  o 


where 


R2  =  r2  +  r  2  -  2rr  cos(e-0  )  +  (z-z  )2 
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